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The functions and elasticities of the cell are largely related to the structures of the cytoskeletons 
underlying the lipid bi-layer. Among various cell types, the Red Blood Cell (RBC) possesses a 
relatively simple cytoskeletal structure. Underneath the membrane, the RBC cytoskeleton takes 
the form of a two dimensional triangular network, consisting of nodes of actins (and other proteins) 
and edges of spectrins. Recent experiments focusing on the malaria infected RBCs (iRBCs) showed 
that there is a correlation between the elongation of spectrins in the cytoskeletal network and the 
stiffening of the iRBCs. Here we rationalize the correlation between these two observations by 
combining the worm-like chain (WLC) model for single spectrins and the Effective Medium Theory 
(EMT) for the network elasticity. We specifically focus on how the disorders in the cytoskeletal 
network affect its macroscopic elasticity. Analytical and numerical solutions from our model reveal 
that the stiffness of the membrane increases with increasing end-to-end distances of spectrins, but 
has a non-monotonic dependence on the variance of the end-to-end distance distributions. These 
predictions are verified quantitively by our AEM and micropipette aspiration measurements of 
iRBCs. The model may, from a molecular level, provide guidelines for future identification of 
new treatment methods for RBC related diseases, such as malaria infection. 


I. INTRODUCTION 

The mechanical properties of a system are largely 
dictated by its structure. The property-structure rela¬ 
tionship has been studied extensively in different fields 
in physics and engineering(e.g., [HE], etc.). In recent 
years, networks in biological systems have drawn much 
attention due to their close relationships to the func¬ 
tions of organisms. Example systems of biopolymer net¬ 
works are cytoskeletons in various cells, which can be 
quasi-1-dimensional (e.g., axons of neuron cells[4l [5]), 
2-dimensional (e.g., red blood cells (RBCs)), or 3- 
dimensional. A biopolymer network can behave very dif¬ 
ferently from a network made of synthesized polymers. 
Firstly, biopolymer networks can be active with the par¬ 
ticipation of ATPs[6]. Secondly, the components of a 
biopolymer network usually follow the worm-like chain 
(WLC) behavior, whose elasticity has an entropic origin 
and thus a nonlinear dependence on the end-to-end dis¬ 
tance of corresponding bio-filaments. 

On the cellular level, the biological functions and 
behaviors of cells are related to their mechanical 
properties[2j El- The mechanical properties of cytoskele¬ 
tons and membranes have been studied intensively in dif¬ 
ferent scenarios, via experiments (e.g., E 0 ID) , simu¬ 
lations (e.g., HDHISl), and theoretical modelings (e.g., 
[3 [H]). In addition, the membrane of RBCs is also 
investigated carefully from both biological and physical 
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perspectives, including the functions of transmembrane 
proteins, and the interactions between protein com¬ 
plexes and the spectrins in the cytoskeletal network (e.g., 
mEsi), etc. As an example to illustrate the relationship 
between the mechanical properties of the membranes and 
the behavior of the cells, earlier studies showed that the 
adhesiveness and hence the mobility of RBCs are strongly 
affected by the stiffness of the membrane (including the 
lipid bilayer and the cytoskeleton) [161 El- Relating to 
this paper, recent experiments reveal that the stiffening 
of RBC after being infected by malaria parasites [IB^ 
HD correlates with the structural transformation in the 
cytoskeleton of the infected RBCs (iRBCs) [22] . A simi¬ 
lar correlation was also observed in iRBCs after chloro- 
quine treatment. In both cases, it is found that, when 
RBCs were infected by malaria or iRBCs were treated 
with chloroquine, the shear modulus of the membrane 
increased with time. Meanwhile, the average length of 
the spectrins that formed the cytoskeletal network was 
increased. It is also noticed that the cytoskeleton mesh 
became more irregular with large holes that were absent 
in the cytoskeleton of normal RBCs, creating a broad dis¬ 
tribution of hole sizes and spectrin lengths. It is known 
that the adhesiveness of RBCs is one of the reasons that 
leads to the fatality of malaria infect ion [23] . Since the 
stiffness of RBCs largely affects the adhesiveness of the 
cells, having a fundamental understanding of the rela¬ 
tionship between the network structure and the macro¬ 
scopic elastic properties will provide us with guidances 
in potentially discovering new drug targets or treatment 
methods. 
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In this study, we focus on rationalizing the correlation 
between the two experimental observations mentioned 
above and investigate how the structural changes at the 
molecular level affect the mechanical properties at the 
cellular level. We present both numerical solutions and 
analytical approximations based on a model combining 
the worm-like chain (WLC) description for single fila¬ 
ments and the Effective Medium Theory (EMT). We use 
the cytoskeletal network of the RBC, approximated by 
a 2-dimensional triangular network, to demonstrate any 
agreement between the model and the experiments. This 
study not only provides a further understanding of the re¬ 
lationship between structures and functions, but also pro¬ 
vides a critical experimental test of the theoretical predic¬ 
tions. In addition, because our model is constructed with 
a general framework, it is expected that it can also be ap¬ 
plied to other scenarios, such as 3-dimensional networks 
of the cytoskeletons in other types of cells, designing of 
new bio-materials, etc. 

II. METHODS 

A. Experiments 

1. Cell culture and enrichment 

The common laboratory P. falciparum 3D7 parasites 
were used for the study. Parasite iRBCs were cultured in 
vitro following a conventional protocol [24j. Culture stage 
synchrony was maintained using sorbitol treatment [25] 
and the early trophozoite stage iRBCs were used for the 
following experiments. 

2. AFM experiment 

Sample preparation: The AEM imaging samples were 
prepared following similar protocol used previously [H]. 
Briefly, synchronized and MACs enriched [25] early 
trophozoite stage malaria iRBCs after different hours of 
1 fiM chloroquine treatment (Ohr, 4hr, 8hr, 16hr and 
24hr) were loaded on PHA-E coated cover slips and then 
incubated for 4 hours to allow sufficient contact time 
for iRBCs to adhere to the substrate and avoid whole 
cell detachment in the shear-wash step. 80 ml of 5P8- 
10 buffer was used to shear-wash the iRBCs adhered to 
PHA-E coated cover slips at an angle of around 20° by sy¬ 
ringe, leaving only the tightly bounded membranes. The 
tight binding between membrane and substrate also max¬ 
imally ensured the structure of the cytoskeleton remained 
in situ. The cytoplasmic-surface-exposed samples were 
checked under phase contrast microscope (Olympus X71) 
and then vacuum dried before imaging using AEM (Eig. 
SI). 

AFM imomng: AEM imaging was carried on JPK 
Nanowizard® AEM with tapping mode (air) using SSS- 
NCHR AEM tip (NANOSENSORSTM) with tip radius 


of 2 nm, and Dimension EastScanTM Atomic Eorce Mi¬ 
croscope (Bruker) using East Scan AEM tip with tip ra¬ 
dius of 5 nm. Images were captured at the resolution 
of 512 X 512 pixels for 10 jam x 10 jam or 1 jam x 1 
fim at scan rate of 0.5 to 1 Hz depending on the scan 
scale and image quality. Around 10 iRBCs were scanned 
for each chloroquine treatment condition and a well con¬ 
nected spectrin area was selected to represent the cell for 
spectrin lengths data collection. 10 to 30 spectrin lengths 
were collected from each cell and averaged to represent 
the spectrin length of that cell. 

Measurement of the spectrin length: The spectrin 
length was measured as described previously [22]. As 
shown in Eig. S2(a) and (b), white lines that were traced 
along the spectrins between two junctions were measured 
as the lengths of the spectrins. 

3. Micropipette experiment 

Different durations of chloroquine treated malaria iR¬ 
BCs were prepared at a ratio of 1: 2000 (cell pellet: 1% 
BSA solution in PBS). 400-600 fil of the diluted sample 
was pipetted into a home-made micropipette aspiration 
microscope cell holder which was mounted on the Olym¬ 
pus X71 microscope stage using masking tape. Borosil- 
icate glass tubings (Sutter) were pulled using a Sutter 
2000 micropipette puller, then forged and cut into mi¬ 
cropipette tip with an inner diameter of 1 - 2 jam using 
Narishige Microforge. A syringe pump was used to apply 
a negative suction pressure to aspirate the iRBC into the 
micropipette. A simple relationship between elongation 
of cell membrane aspirated into the micropipette and the 
applied negative pressure is given by HE]: 

(AP X Dp)/ia = 2.45 x ALp/Dp (1) 

with ALp/Dp > 1, and AP is the applied negative pres¬ 
sure, Dp is the micropipette inner diameter, /r is the 
membrane shear elastic modulus and ALp is the aspi¬ 
rated length of the cell membrane into the micropipette. 
Series of images were captured using a CCD camera ev¬ 
ery second. By measuring the aspirated length and the 
applied negative pressure, the membrane shear modulus 
was obtained. 


B. Irregular network of worm-like chains (WLC) 

Normal RBC possesses a triangular cytoskeletal net¬ 
work consisting of nodes of protein complexes and edges 
of spectrins. In our theoretical study of the RBC cy¬ 
toskeleton, we regard the network formed by equilateral 
triangles as the regular network. We treat this regular 
network as our reference state and study how disorders 
added to the network affect its elasticity. 

Recent experiments [22] found that in malaria infected 
RBCs (iRBCs) or iRBCs after drug (chloroquine) treat¬ 
ment, two types of disorders are introduced to the RBC 
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Figure 1. Sketches of the model, (a) A sketch of the 2D tri¬ 
angular network with disorders, resembling the cytoskeleton 
structure of the RBC (or iRBC). The lengths Zij between two 
nodes are non-uniform and some connections are missing com¬ 
paring to a regular network, (b) Based on the WLC model 
(Eq. and Eq. |^, the length variation in zij is mapped to 
the variation in the equivalent spring constant kij for small 
perturbations, (c) An EMT is applied to reduce the disor¬ 
dered network to a regular network with uniform spring con¬ 
stant km , from which the macroscopic in-plane shear modulus 
can be calculated. 


introduces a topological disorder to the network. Here 
we use a parameter p to represent the probability with 
which an edge of spectrin is present in the regular net¬ 
work (then 1—p represents the probability with which an 
edge of spectrin is missing from the regular network). For 
example, p = 1 corresponds to the complete network, and 
p = 2l3 corresponds to an average out-degree of 4. In our 
model, we generally treat p as an independent variable. 
But it is possible that p can be a function of stresses in 
the network as suggested by previous studies focusing on 
the binding/unbinding kinetics between ligands and re¬ 
ceptors (e.g., [26]). This case is further discussed in the 
Discussion section. 

For individual spectrins, it is well-known that the elas¬ 
ticity has an entropic origin. Here, we apply the worm¬ 
like chain (WLC) model that is used to study semi- 
flexible polymers (e.g., [27]-[29|) to single spectrins, which 
basically gives us a nonlinear force-extension curve [77]. 
Comparing to the end-to-end distance of the spectrins 
in the network (^ 50 nm), the persistence length {Ip ^ 
lOnm) of the spectrin is very short. Hence we need not 
consider the effect of the spectrins’ bending rigidity on 
the network. For a brief review, in the WLC model, the 
persistence length Ip describes the length scale, below 
which the bending energy dominates the thermal excita¬ 
tions. From the WLC model, the single spectrin follows 
the force-extension relation as shown below [77]: 



1 1 

4(1 -z/L)2 “ 4 


(3) 


cytoskeletal network (Fig. [^. Firstly, the average end- 
to-end distance of the spectrins increases, which is di¬ 
rectly measured in our experiment (Fig. |^(b)). Secondly, 
the variance of the end-to-end distance distribution is 
also much larger than that in the regular triangular net¬ 
work (Fig. |^(b) and [22]), probably due to the absence 
of some spectrin links after infection or drug-treatment. 
To perform a quantitative investigation of how these two 
changes affect the stiffness of the RBC network, we in¬ 
troduce the disorders in terms of a Gaussian distribution 
of the spectrins’ end-to-end lengths 2 : in the network, for 
which the probability density function can be written as: 

V{z) = ^ e 2 <t 2 0 < z < L (2) 

V27rcr 

Here, L is the natural length of the spectrins following the 
WLC description elaborated later, zq is the average end- 
to-end length of the spectrins, is the variance of the 
end-to-end length distribution, and (7 is a normalization 
constant. 

Each node in the network formed by protein complexes 
defines a point connecting the cytoskeleton and the mem¬ 
brane. In the regular network, each node is linked by 6 
spectrins, i.e., an out-degree of 6. But in reality, some 
spectrins may lose the connections with the nodes, which 


where p = IjksT is the Boltzmann factor with ks the 
Boltzmann constant and T is the temperature, z is the 
end-to-end distance of the spectrin and L ^ 200 nm is 
the natural length of the spectrin that is treated as in- 
extensible. With L ^ 200 nm and Ip ^ 10 nm^ the 
root mean square end-to-end distance of a spectrin free 
of stress is about 60 nm, comparable to observed values. 
In general, the spectrins in the RBC cytoskeleton has an 
average end-to-end distance ranging from about 50 nm 
to 100 nm. With small perturbations to the spectrin 
from its equilibrium position, we can expand the force- 
extension relation (eq. (§) to first order and obtain an 
effective spring constant k for the spectrin, which is a 
function of the end-to-end distance 2 : of the spectrin at 
equilibrium. Simply, this extension-related spring con¬ 
stant k is obtained by taking the derivative of the force 
with respect to the extension, which gives us 

fc(^) = g = (i + i(i-|)-3)//?;,L (4) 

Importantly, we should note here that, even the effec¬ 
tive spring constant k depends on the end-to-end distance 
2 : nonlinearly, its effect in our model is simply to setup 
an initial equilibrium distribution of the spring constants 
across the entire network according to the distribution 
of the end-to-end distance of the spectrins. Under this 
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assumption, we focus on the linear response of the disor¬ 
dered network (small perturbations), in which case each 
spectrin is treated as a Hookean spring with a distinct 
spring constant and then we apply an Effective Medium 
Theory (EMT) to evaluate the elasticity of the disordered 
network. 


C. Effective medium theory 

To understand how the disorders in the RBC cytoskele- 
ton affect its stiffness, we adopt an Effective Medium 
Theory (EMT) approach here. In previous studies, EMT 
is used to understand the elastic properties of different 
types of disordered networks. Examples include but are 
not limited to central-force spring networks [301 EH, 2D 
networks under tension [32H34] . networks at large defor¬ 
mations [35], networks formed by filaments with finite 
bending rigidities IM1E2, networks with nonlinear cross¬ 
linkers [38], etc. Basically, the EMT maps the disordered 
network to an equivalent regular network and extract 
elastic constants from the regular network. To construct 
the equivalent regular network, the springs in the network 
with distinct spring constants are replaced by springs 
with the same spring constant km self-consistently. The 
self-consistency requires that the average extra displace¬ 
ment Su caused by this replacement procedure must be 
0 over the entire network (< 6u >= 0) [30]. It leads to 
the following equation for the effective spring constant in 
the regular network [30]: 




without treatment 24 hrs after treatment 


Figure 2. Experimental results, (a): Shear modulus fi as 
a function of the time after chloroquine treatment, (b): The 
end-to-end distance z as a function of the time after the treat¬ 
ment. In both (a) and (b), 0 hour data corresponds to cells 
without chloroquine treatment. The shaded areas between 
the blue and orange curves in (a) and (b) indicate the 90% 
point-wise confidence band for the mean values respectively, 
(c) and (d): Typical images from AFM at different stages of 
the treatment (Without chloroquine treatment for (c) and 24 
hrs after the treatment for (d)). The color indicates the rel¬ 
ative height of the corresponding point. The scale bar in (c) 
and (d) is 1 /xm. 


where V{k) describes the distribution of the spring con¬ 
stants in the network, a is related to the topology of and 
pre-stresses in the network [32] and can be understood as 
follows: When two adjacent nodes in the regular network 
is displaced with respect to each other, the response of 
the network can be described by an effective spring con¬ 
stant K' = kmjoL^ where 0 < a < 1 takes into account 
the contribution of the whole network. Eq ^ can be 
solved analytically for specific probability density func¬ 
tions V{k) and numerically for other forms of V{k). 

III. RESULTS 

A. Experimental measurements 

In the cytoskeleton of normal RBCs, the average end- 
to-end distance of the spectrins, i.e., the average distance 
between nodes, is about 50 — 60 nm. Here we summarize 
our experimental measurements using Atomic Eorce Mi¬ 
croscopy (AEM) and micropipette aspiration techniques. 
The measurements show correlation between the stiffen¬ 
ing of the iRBC after chloroquine treatment as well as 


elongation of spectrins which motivate our theoretical 
model. 

The shear modulus of the iRBC is measured by mi¬ 
cropipette aspiration [3^ |4Q] . Eig. [^a) shows that the 
shear modulus of iRBC increases by about a factor of 2 af¬ 
ter 24 hours of chloroquine treatment. The cytoskeleton 
structure is measured using AEM with a similar method 
as that in [22]. Erom the AEM images, we measured 
the end-to-end distances of spectrins in the cytoskeleton 
network (Eig. (b) and Eig. S3). The result shows an 
overall increase of the end-to-end distances after chloro¬ 
quine treatment. The typical AEM images are shown in 
Eig. I^c) and (d). The larger dark spots in Eig. i(d) 
indicate larger holes in the cytoskeletal network, consis¬ 
tent with our measurement of the end-to-end distances 
of spectrins. 

B. Numerical calculation 

Eor a general probability density function V{k)^ Eq. 
can be solved numerically. Here we consider the varia¬ 
tion of the lengths of spectrins in the network, as well 
as the probability of missing spectrins between nodes. 
Because the effective spring constant of the spectrin de¬ 
pends on the end-to-end distance, the length distribution 
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introduced by Eq. induces the probability distribution 
of elastic constants for spectrins. Combined with the 
probability of missing a spectrin connection (1 — p), the 
probability density function of the spring constants in the 
cytoskeletal network can be written as 

~ d z 

■p(^k)=p*V{z)— + {l-p)S{k) (6) 

where V{z) is the distribution of the end-to-end distance 
^ and normalized for 0 < 2 : < L, i.e., P(z)dz = 1, 
and S(k) is the Dirac Delta function. In the first term on 
the right-hand-side, V{z)j^ gives the distribution of the 
spring constant k. Using the probability density function 
(eq. i the EMT equation (Eq. is re-written as 

As discussed earlier, the value of a depends on the geom¬ 
etry of and the pre-stresses in the network. Here we take 
the value a = 2/3 for stress-free triangular network [30] 
and solve km numerically. Tang et.al. [32| shows the rela¬ 
tion between a and pre-strain in the network. We found 
that changing the value of a (from 1/3 to 2/3) changes 
the value of km but does not change the qualitative de¬ 
pendence of km on p, zo, and a as shown in Eig. Once 
we obtain the effective spring constant km for the equiv¬ 
alent regular network, the in-plane shear modulus to 
the lowest order, is linearly proportional to km- Hence, 
to study how p depends on the disorders, it is equivalent 
to find out the dependence of km on the connection for¬ 
mation probability p, the average end-to-end distance zq, 
and the variance o^ z. 

Our numerical results are summarized in Eig. With 
fixed value of p = 0.8 (Eig. |^a)), the effective spring 
constant km increases with a for relatively small values 
of 2 ^ 0 , but decreases with a for relatively large values of zq 
{zq > 100 nm). On the other hand, km always increases 
with increasing zq given the value of a. This supports the 
prediction that the experimentally observed cell stiffening 
is related to the observed lengthening of spectrins in the 
cytoskeleton network. With fixed value of a = 30 nm 
(Fig. ib)), as one may expect, km increases with in¬ 
creasing Zq given the value of p, and also increases with 
increasing p given the value of zq. With fixed value of 
Zq = 100 nm (Eig. I^c)), km increases with increasing 
p given the value of a. But similar as shown in Eig. 
ia), km increases with increasing a for relatively high 
values of p (p > 0.86) while it decreases with increasing 
a for relatively low values of p. Therefore, our numeri¬ 
cal results suggest that, given the value of p (e.g., Eig. 
[^a)), for relatively “strong” network (large values of km)^ 
increasing the variance of 2 : tends to weaken the net¬ 
work, while for relatively “weak” network (small value of 
km), increasing the variance tends to stiffen the net¬ 
work. However, given a relatively medium value of zq 
(e.g., Zq = 100 nm in Eig. |^c)), increasing a stiffens rel¬ 
atively “strong” network, but weakens relatively “weak” 
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Figure 3. Numerical results for the effective spring constant 
km- (a) The dependence of km on zo and a (Eq. with the 
probability for one spectrin edge to present p = 0.8 fixed, (b) 
The dependence of km on zq and p with cr = 30 nm fixed, (c) 
The dependence of km on a and p with zq = 100 nm fixed. 
The color bars show the value of km in the unit of fiN/m. 
The units match the units used in relevant experiments. 

network. This nontrivial dependence of km on a may be 
utilized in cellular functions or material designs. 


C. Analytical solutions 

Here we adopt analytical solutions from previous stud¬ 
ies (e.g., [3Ql 137] ) but replace the harmonic springs with 
worm-like chains. To make the model traceable, we only 
consider the deviation from the regular triangular net¬ 
work caused by missing spectrin connections and neglect 
the length heterogeneity. Aiming at the dependence of 
shear modulus on the average spectrin length, here we 
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will establish a correlation between the experimentally 
observed spectrin lengthening and cell stiffening phenom¬ 
ena. In this case, the distribution function takes the fol¬ 
lowing form [30l |35] 

V{k')=p6{k' - k{z)) + {l-p)5{k') 

Different from purely linear springs, stretching spectrins 
not only introduces internal pre-stresses to the network, 
but also increases the effective spring constants k{z) of 
individual spectrins. If we use the average extension < 
z > to replace the inhomogeneous end-to-end distance 
Eq. can be solved exactly [30]. Again, the length- 
dependent spring constant only sets up the initial value 
of k{z). After that, we treat ah the spectrins as Hookean 
springs and k as independent of z, and only consider the 
linear response of the network for small perturbations. 
With this consideration, Eq. [^is solved to give us: 


km — 


P-Pc 


^-Pc plpL 






Pc 


( 8 ) 


Here we consider p as a constant independent of initial 
stresses. The case of a stress-dependent p will be dis¬ 
cussed in the Discussion section. Thus, we map the disor¬ 
dered network to a regular triangular network consisting 
of springs with a single spring constant km- It is clear 
that, in our model, cell stiffening comes directly from 
the nonlinearity of the WLC model, as shown by the z- 
dependent term on the right hand side of Eq. The 
macroscopic in-plane shear modulus p of the 2D network 
is linearly proportional to the effective spring constant, 
i.e., p (X km- 

On the other hand, it is known that pre-stresses can 
stiffen the network (e.g., [sa El [35]), a similar effect as 
discussed in subsection Numerical calculation. In a cel¬ 
lular environment, the pre-stresses can be caused by var¬ 
ious processes, e.g., the expression and export of trans¬ 
membrane proteins onto the iRBC membrane that bind 
to spectrins, the formation of “knobs” at dispersed sites 
on the RBC membrane that can link to near-by spec¬ 
trins, the expansion of the cell’s volume due to the in¬ 
vaded and multiplying malarial parasites and so on. In 
our model, without knowing more details about these 
intracellular processes, we use a single stress-dependent 
pre-factor a{P) to account for the influence of the pre¬ 
stresses on the elastic properties of the cytoskeletal net¬ 
work. With that in mind, the shear modulus of the 2D 
network is given by 


P — (^{P^km 


<P) 

PlpL 


(1 + ^(1 


Z.-3) P-Pc 

P P-Pc 


(9) 


where P is the pre-stress in the network, defined to be 
negative for networks under tension. 

Fig. i a) shows our analytical results for different val¬ 
ues of Pc- We focus on the dependence of the shear mod¬ 
ulus p on the average end-to-end distance zq and we use 



Figure 4. Analytical results and comparison with experi¬ 
ments. (a) Dependence of the shear modulus p on the av¬ 
erage end-to-end distance zq. Red dots with error bars are 
experimental data. Different curves correspond to different 
values of pc- Inset: comparison between the analytical re¬ 
sult and the experiments. The black curve is the analytical 
result, taking the same parameters as the black curve in (a) 
{Pc = 0.45). The shaded area indicates the 90% point-wise 
conhdence band for the experimentally obtained mean values 
of the shear modulus, (b) Non-m onot onic dependence of p 
on zo for different values of ^ (Eq. when the unbinding 
kinetics is taken into account. 


p = 0.55 that is approximated from our AEM measure¬ 
ments. Qualitatively, it is clear that the shear modulus of 
the network increases with increasing average end-to-end 
length Zq of the spectrins in the network. By adjusting 
the value of Pc and a(P), our results agree with experi¬ 
ments (red dots with error bars in Eig. [^a)) well. The 
comparison with experiments is elaborated below. 


D. Comparison with experiments 

Eq. is used to compare with experiments. By us¬ 
ing micropipette aspiration techniques, we found that the 
malaria iRBCs become stiffer after chloroquine treatment 
(Fig. i (a)). In the meantime, through AEM measure¬ 
ments, we found that the average end-to-end distance 
of spectrins in the cytoskeleton increases after the drug 
treatment (Eig. (b)). This coincidence is also no¬ 
ticed in iRBCs at different stages of malaria infection 
[TBHll]. In our experiments, both the shear modulus p 
and the spectrin lengths were measured (Eig. |^a)), and 
the data showed a positive correlation between these two 
quantities. In the meantime, the average connectivity is 
also obtained, which is around 3.3 across several sam¬ 
ples. So we use this number to determine the value of 
p ~ 3.3/6 = 0.55 in our model. Eor other parameters, we 
use L = 200 nm, T = 300 K and Ip = 10 nm, which agree 
with the values reported in the experiments or literatures 
( [411143] ). a and Pc are used as fitting parameters. Using 
the values a = 15 andpc = 0.45, our results (Eq.agree 
well with the experimental observations (Eig. ^a) and 
the inset). Eor the black curve in the inset, the mean 
values of the lengths of spectrins at different times are 
used in Eq. [^to obtain the values of p respectively. Due 
to the relatively large error bars in the experiments, the 
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values of a and Pc can vary in a certain range. Better 
estimations of the parameter values are worthy of fur¬ 
ther investigations and more precise experimental mea¬ 
surements in our future work. However, the agreement 
supports the prediction that the shear modulus of the 
network is largely affected by the average length of edges 
(e.g., spectrins) in the cytoskeletal network. 


IV. DISCUSSION 


Several questions deserve further discussions here. 
Firstly, the binding/unbinding kinetics between spec¬ 
trins and the nodes (protein complexes) can be force- 
dependent. Here we simply assume an exponential de¬ 
pendence of the edge missing probability on the force, 
i.e., (1 — p) ^ exp(;dx^/), where / = f(z) depends on 
the end-to-end distances of the spectrins as in Eq. H By 
plugging in p = 1 — ^exp(/3a;’^/(2)) in Eq. we obtain: 


= 


“(-P)n , -Cexp(^a;’^/(2:))) -Pc 


f3lpL 




Pc 


( 10 ) 

where is a parameter describing the length scale be¬ 
tween the binding potential minimum and maximum. 
In Eq. 2 ; is approximated as the average end-to- 

end distance < z >. Two competing effects are im¬ 
mediately noticed from Eq. 10 The (1 + |(1 — 
term coming from the WLC model predicts cell stiffen¬ 
ing due to lengthening of spectrins, while the other term 
1 — f{z)) weakens the network. Eor relatively 

short spectrin lengths, the cell stiffens due to the stiffen- 
ings of spectrins as their lengths are increased. However, 
at large spectrin lengths that may be beyond the mea¬ 
surements of current experiments, the weakening term 
dictates the elastic properties of the network and the 
bond forming probability p may even drop below the 
rigidity percolation threshold Pc- Thus a turning point 
should be observed if the length of the spectrins keep in¬ 
creasing, dividing the two different dependences of shear 
modulus on 2 ;. Eig. [^b) confirms this prediction, and 
shows the dependence of the shear modulus p on the av¬ 
erage end-to-end distance zq of the spectrins for different 
values of In Eig. [^b), we use x^ = 1 nm that is in the 
range of values for similar problems [ 16 ] and Pc = 0 . 45 . 
Other parameters take the same values as those for Eig. 
I^a). The curves in Eig. |^b) are slightly shifted from 
those in Eig. [^a) due to the differences between Eq. 
and Eq. The agreement with experiments (red dots 
with error bars) can be retained by only slightly chang¬ 
ing the value of a (a = 12 for the gray dashed curve in 
Eig. I^b)) with other parameters stay the same. Having 
that said, here we are more interested in the qualitative 
non-monotonic dependence of p on zq^ while the exact 
shape of the curve needs to be further determined by 
parameters from future experimental investigations. 

Secondly, it is noticed that pre-stress can significantly 
alter the geometry or topology of the network through 


stretching of spectrins. The effects of pre-stresses in cells 
and cytoskeletons have been studied extensively previ¬ 
ously regarding the stability of the triangular cytoskele¬ 
tal network and the origin of the stresses are worthy of 
further studies. One scenario is the stress generated in 
the cytoskeleton of malaria infected RBCs. It is known 
that after the malaria infection, “knobs” are expressed 
across the membrane, which are complexes of malarial 
parasite exported proteins. These knobs are found ca¬ 
pable of binding to spectrins and can possibly introduce 
stresses in the cytoskeletal network[22]. A different sce¬ 
nario is the iRBCs after chloroquine treatment. In this 
case, the oxidative stress resulted from the drug treat¬ 
ment may be responsible for the change of mechanical 
properties of the cytoskeleton [9]. We suspect that this 
oxidative stress may be related to the formation of new 
protein complexes that binds to spectrins, and subse¬ 
quently generate stresses in the network. In this case, 
the effect of chloroquine treatment is partially impaired 
by the continuous stiffening of the cells, because stiffer 
cells are sometimes easier to block the blood circulation. 
Eurther experimental investigations of the origin of the 
stresses are our main focus in the future and will improve 
our understanding of the mechanism of fatal diseases and 
drug treatments. 

Erom a broader perspective, although here we looked 
at the 2-dimensional network, especially the RBC cy¬ 
toskeleton as an example, our model in general bridges 
the properties of biopolymers at the molecular level and 
the elastic properties at the cellular level. Similar models 
may be extended to other dimensions and shapes. Our 
results also show a possible mechanism for cells to adjust 
their stiffness from the molecular level and may stimulate 
new ideas in material designs. 


V. CONCLUSION 


To conclude, we investigated the relationship between 
the network structure and its stiffness. We focused on 
how the disorders introduced to a regular 2D triangu¬ 
lar network affect its macroscopic in-plane shear mod¬ 
ulus. By applying the worm-like chain model to single 
spectrins that form the edges of the network, we found 
that the shear modulus of the network increases with the 
average length of the spectrins. Our result agrees with 
previous studies in that removing some edges with prob¬ 
ability 1 — p weakens the network, which shows a com¬ 
peting effect with the stiffening caused by lengthening 
of spectrins (worm-like chains). On the other hand, as 
discussed in the results section, depending on the val¬ 
ues of the average length zq or p, the shear modulus 
may increase or decrease with increasing variance (cr^) 
of the length distribution (Eig. [^a)(c)). Eurthermore, if 
the stress-dependent binding/unbinding kinetics is taken 
into account, the shear modulus has a first-increase-then- 
decrease behavior as the average length of the spectrins 
increases. More importantly, our results agree well with 





recent experimental observations for malaria infected red paper, 
blood cells (iRBCs), which shows a correlation between 
the lengthening of spectrins (edges of the cytoskeleton 
network) and stiffening of the cells (increase of shear 
modulus). Further investigations of what causes the 
lengthening of spectrins or the internal stresses in the cy¬ 
toskeleton will be valuable in identifying new treatment 
methods or targeted therapy for RBC-related diseases, 
such as malaria infection. 

ACKNOWLEDGMENTS 


AUTHOR CONTRIBUTIONS 

We acknowledge the financial assistance of Singapore- 
LL, XX, CTL, and JC designed the research. LL and MIT Alliance for Research and Technology (SMART) 

XX performed the research. LL, CTL, and JC wrote the and National Science Foundation (NSF) (CHE-1112825). 


[1] D. Discher, N. Mohandas, and E. Evans, Science 266, 
1032 (1994). 

[2] S. Suresh, Journal of Materials Research 21, 1871 (2006). 

[3] C. P. Broedersz and F. C. MacKintosh, Rev. Mod. Phys. 
86, 995 (2014). 

[4] K. Xu, G. Zhong, and X. Zhuang, Science 339, 452 

( 2012 ). 

[5] L. Lai and J. Cao, The Journal of Chemical Physics 141, 
015101 (2014). 

[6] J. H. Kim and J. Cao, in preparation. 

[7] S. Huang, A. Amaladoss, M. Liu, H. Chen, R. Zhang, 
P. R. Preiser, M. Dao, and J. Han, Infection and Immu¬ 
nity 82, 2532 (2014). 

[8] G. Lenormand, S. Henon, A. Richer!, J. Simeon, and 
F. Gallet, Biophys J. 81(1), 43 (2001). 

[9] J. P. Hale, C. P. Winlove, and P. G. Petrov, Biophys J. 
101, 1921 (2011). 

[10] M. J. Saxton, Biophys J. 57, 1167 (1990). 

[11] D. H. Boal, Biophys J. 67, 521 (1994). 

[12] J. C. Hansen, R. Skalak, S. Chien, and A. Hoger, Biophys 
J. 72(5), 2369 (1997). 

[13] S. K. Boey, D. H. Boal, and D. E. Discher, Biophys J. 
75, 1573 (1998). 

[14] N. S. Gov, New Journal of Physics 9, 429 (2007). 

[15] T. J. Byers and D. Branton, Proceedings of the National 
Academy of Sciences 82, 6153 (1985). 

[16] A. Efremov and J. Cao, Biophys J. 101, 1032 (2011). 

[17] X. Xu, A. K. Efremov, A. Li, L. Lai, M. Dao, C. T. Lim, 
and J. Cao, PLoS ONE 8, e64763 (2013). 

[18] S. Chien, K.-L. P. Sung, R. Skalak, and U. Shunichi, 
Biophys. J. 24, 463 (1978). 

[19] G. Nash, E. O’Brien, E. Cordon-Smith, and J. Dor- 
mandy. Blood 74(2), 855 (1989). 

[20] M. Paulitschke and G. B. Nash, J Lab Clin Med. 122, 
581 (1993). 

[21] F. K. Glenister, R. L. Coppel, A. F. Cowman, N. Mo¬ 
handas, and B. M. Cooke, Blood 99, 1060 (2002). 

[22] H. Shi, Z. Liu, A. Li, J. Yin, A. G. L. Chong, K. S. W. 
Tan, Y. Zhang, and C. T. Lim, PLoS ONE 8, e61170 
(2013). 

[23] K. R. Hughes, G. A. Biagini, and A. G. Craig, Molecular 
& Biochemical Parasitology 169, 71 (2010). 

[24] W. Trager and J. Jensen, Science 193, 673 (1976). 


[25] C. Ribaut, A. Berry, S. Chevalley, K. Reybier, 1. Morlais, 
D. Parzy, F. Nepveu, F. Benoit-Vical, and A. Valentin, 
Malaria Journal 7, 45 (2008). 

[26] G. Bell, Science 200, 618 (1978). 

[27] J. F. Marko and E. D. Siggia, Macromolecules 28, 8759 
(1995). 

[28] S. Yang, J. B. Witkoskie, and J. Cao, Chemical Physics 
Letters 377, 399 (2003). 

[29] X. Xu, B. J. R. Thio, and J. Cao, The Journal of Physical 
Chemistry Letters 5, 2868 (2014). 

[30] S. Feng, M. Thorpe, and E. Garboczi, Phys. Rev. B 31, 
276 (1985). 

[31] M. Thorpe and E. Garboczi, Phys. Rev. B 42, 8405 
(1990). 

[32] W. Tang and M. Thorpe, Phys. Rev. B 37, 5539 (1988). 
33 D. Boal, U. Seifert, and J. Shillcock, Phys. Rev. E 48, 

4274 (1993). 

[34] D. H. Boal, Biol. Bull. 194, 331 (1998). 

[35] M. Sheinman, C. Broedersz, and F. MacKintosh, Phys. 
Rev. E 85, 021801 (2012). 

[36] M. Das, F. MacKintosh, and A. Levine, Phys. Rev. Lett. 
99, 038101 (2007). 

[37] X. Mao, O. Stenull, and T. Lubensky, Phys. Rev. E 87, 
042601 (2013). 

[38] C. Broedersz, C. Storm, and F. MacKintosh, Phys. Rev. 
E 79, 061914 (2009). 

[39] E. Evans, Biophysical Journal 13, 941 (1973). 

[40] R. M. Hochmuth, Journal of Biomechanics 33, 15 (2000). 

[41] B. T. Stokke, A. Mikkelsen, and A. Elgsaeter, Biochim- 
ica et Biophysica Acta 816, 102 (1985). 

[42] K. Svoboda, C. F. Schmidt, D. Branton, and S. M. 
Block, Biophys. J. 63, 784 (1992). 

[43] J. Li, M. Dao, C. T. Lim, and S. Suresh, Biophys. J. 88, 
3707 (2005). 



1 


SUPPORTING FIGURES 


/a 


MACS enriched 
iRBCs 



^^P ^ PhA-E 

Cover slip 

V_ J 







Attached 



membrane 

/ 


1“ 





J 



Figure SI. Cytoplasmic-surface-exposed membrane preparation process. A. MACs enriched iRBCs was first adhered to a 
PHA-E coated cover slip. B. Shear washing of adhered iRBCs using the 5P8-10 buffer at an angle of around 20° by syringe to 
peel off the top iRBC membrane. C. Cytoplasmic-surface-exposed membrane remained on the cover slip after shear washing. 
D. Phase contrast imaging of iRBC cytoplasmic-surface-exposed-samples. 



Figure S2. (a): 3D structure spectrin network constructed by JPKPM Data processing software for AFM image obtained 

using JPK Nanowizard(R)AFM. (b): 3D structure spectrin network constructed by using NanoScope Analysis 1.5 software for 
AFM images obtained using Dimension FastScanTM Atomic Force Microscope. 
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Figure S3. Distribution of spectrin lengths measured after different durations of chloroquine treatment. The time durations 
of treatment are shown in the legend. 













